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Abstract. A possible signal of new phenomena emerging in the global characteristics 
of multiparticle production in hadron interactions at TeV energies is studied. The 
multiplicity distributions of charged particles measured in proton-proton collisions at 
the Large Hadron Collider at CERN are analyzed in the phenomenological framework 
of the weighted superposition of three negative binomial distribution functions. The 
examination of the experimental data indicates the existence of a narrow peak at low 
multiplicities which can be described by a separate component in the total distribution. 
The multiplicity characteristics of the third component reveal approximate energy and 
pseudorapidity invariance, of which a physical explanation represents a challenging 
problem in high energy multiple particle production. 
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1. Introduction 

The multiple production of particles in high energy hadron interactions has received 
a lot of attention over the years. A renewal in interest in this topic coincided with 
the operation of the Large Hadron Collider (LHC) at CERN. Among the first results 
obtained at the LHC were the multiplicity measurements in proton-proton collisions. 
The ATLAS, CMS, ALICE and LHCb Collaborations provided high statistic data on 
multiplicity distributions (MDs) of charged particles produced in the new super high 
energy domain. 

The study of particle production as a function of multiplicity has revealed the very 
popular Koba-Nielsen-Olesen (KNO) scaling [1]. In pp/pp collisions, scaling in the full 
phase-space holds up to the highest energy of the CERN Intersecting Storage Rings but 
it is clearly violated |2] in the energy region of the CERN Super Proton Synchrotron 
collider and beyond. The strong violation of KNO scaling was observed also at the 
LHC energy y/s = 7 TeV in the limited pseudorapidity intervals {\ri\ < 2.4) though 
for small pseudorapidity windows {\ri\ < 0.5) the scaling is approximately valid [3J. A 
related phenomenon is the so-called negative binomial regularity which is the occurrence 
of the negative binomial distribution (NBD) in different interactions over a wide range 
of collision energies. The UA5 Collaboration showed that MDs of charged particles 
in non-single-diffractive (NSD) pp/pp collisions can be described by NBD up to the 
energy = 546 GeV both in the full phase-space |1] and in symmetric pseudorapidity 
windows [5]. First analysis of MDs measured by the ALICE Collaboration [6] indicated 
that NBD describes the data in the small pseudorapidity window \ri\ < 0.5 up to the 
energy y/s = 2360 GeV [7]. The accuracy of such a description is however subject 
to further disputes and investigations [8]. Much effort has been made to explain the 
negative binomial form of MDs observed in many situations, still its physical origin has 
not been fully understood [5]. 

The shape of MD of charged particles produced in sufficiently large windows of 
pseudorapidity in high energy hadron-hadron collisions is quite different. The full 
phase-space data on multiplicities [10] obtained from the NSD events in pp collisions 
at a/s = 900 GeV indicated that besides KNO, the negative binomial regularity is 
violated as well. The measurements of MD at the energy ^/s = 1800 GeV by the E735 
Collaboration [11] at the Fermilab Tevatron showed even stronger deviation of data 
from NBD. Though with large mutual discrepancies at high multiplicities, both data 
demonstrate a narrow maximum and a shoulder-like structure around n ~ 2 < n >. 
Despite its correct qualitative behavior, single NBD is not sufficient to describe the 
experimental data. More functions were suggested [121 [13] to investigate the properties 
of high energy multiple particle production. A systematic study of the complex form of 
MDs was performed in the framework of the two-component model [H] and its three- 
component modification [15]. The multicomponent quark gluon string model (QGSM) 
[T6] was used [17] to study the role of multiparton coUisions. There also exist other 
approaches to particle production at high energies (for a review see [9]). 
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The ATLAS [IB], CMS [5] and ALICE [IS] Collaborations measured MDs of charged 
particles in proton-proton collisions at the LHC energies a/s = 0.9, 2.36 and 7 TeV 
using different kinematic ranges and different classes of events. The measurements in 
the central pseudorapidity windows allow to study evolution of the narrow maximum 
and the broad shoulder of the distributions at high n both with collision energy and 
pseudorapidity. The shoulder-like structure of MD was investigated in the framework 
of independent pair parton interaction (IPPl) [20j . The IPPI model was applied pTj to 
experimental data measured by the CMS Collaboration [3J at y/s = 7 TeV in various 
pseudorapidity intervals. The CMS data were analyzed [22] by weighted superposition 
of two NBDs representing soft and semihard components of particle production. 

In this paper we study the three-component description of data on MD of 
charged particles produced in proton-proton collisions at the LHC. Each component 
is represented by a single NBD and parametrized by two parameters n and k. The total 
distribution is the weighted sum of three NBDs with eight free parameters altogether. 
The high statistic data on MD measured at the LHC allow us to extract these parameters 
in most cases reasonably well. The analysis of the ATLAS measurements [18| shows 
us that the two-component description of MD is unsatisfactory. Within the multi- 
component NBD parametrization, the structure of the data indicates the necessity of a 
third component in the region of low multiplicities. The shape of the ATLAS data at low 
n is critical to the IPPI model and to the two-component superposition of NBDs. The 
goal of our analysis is to obtain a detailed description of the high statistic data on MD 
in pp collisions at the LHC in the framework of weighted superposition of three NBDs 
and study its properties. We demonstrate that, besides the shoulder-like structure at 
high n, the data manifest a distinct peak at a maximum which can be described by the 
third component of the total distribution. We show that the obtained parametrization 
of the LHC data on MD reveals some invariant features. 



2. Weighted superposition of NBDs 



The multiple production of particles in the soft processes at high energies includes 
highly non-perturbative effects of QCD. A comparison of experimental data with various 
distributions relies mostly on phenomenological approaches in this field. The prominent 
role in the phenomenological description of data on MD in pp/pp collisions is the two- 
parameter NBD 
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The parameter n is the average multiplicity and k characterizes the width of the 
distribution. The Poisson distribution is obtained for k = oo. A remarkable property 
is the convolution of NBDs which results again in NBD with modified parameters. A 
superposition of NBDs is a well-known idea exploited by different approaches to multiple 
particle production in hadron collisions. It is based on the assumption that MD of the 
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created particles can be expressed as 

N N 



(2) 




where P{n, fii, ki) is given by ([T]). In the general case, the corresponding number of free 
parameters is 3A^ — 1. In specific models of this type, the number of parameters can 
be smaller. In the IPPI model there are three parameters, a, m and k, given by the 
conditions = a*, fii = im and ki = ik. The parameter a is defined by the number 
of the independent active parton pairs via the normalization of the probabilities a, 
to unity. The value of = 6 was used in the analysis [21] of the CMS data [3j at the 
energy a/s = 7 TeV. In the QGSM model the probabihties a, are proportional to the 
cross sections for simultaneous production of 2i chains created in i-fold parton collisions 
[17] . The MD of the i cut-pomerons contribution, P{n,ni, ki), is assumed to have a 
Poisson-like form. The characteristic feature of both models is an increase of the mean 
multiplicities Ui = mi {fii = i[a + 61n(s/soi^)]) with decreasing values of a^. In this 
paper we consider the function ([2]) for = 3 in a formal way and try to establish values 
of the eight independent parameters {fii, ki, a2 and as) to obtain a detailed description 
of data on MD measured with high accuracy in the experiments at the LHC 

3. Analysis of data 

The MDs of charged particles produced in proton-proton collisions at the LHC were 
measured both in the central and forward pseudorapidity region. The phenomenological 
studies of the measurements have been mainly carried out [211 EZl |23l El] with the CMS 
data. In this paper we present results of a combined analysis of data on MD obtained 
by the ATLAS, CMS and ALICE Collaborations in the central interaction region and 
by the LHCb Collaboration in the forward region. 

The ATLAS Collaboration measured [18] the charged particle MD in different 
phase-space regions using various multiplicity cuts in the pseudorapidity window |?7| < 
2.5. The distributions are based on the experimental analysis of a large number of 
events which accumulated in particular at the energies ^/s = 7 and 0.9 TeV. The most 
inclusive phase-space region covered by the measurements corresponds to the conditions 
Pt > 100 MeV/c and rich > 2. The record number of the analyzed events in this region 
exceeds lOM and 0.3M at ^/s = 7 and 0.9 TeV, respectively. A similarly large number 
of events was considered in a data sample with the higher transverse momentum cut 
Pt > 500 MeV/c and rich > 1- The measurements in both regions give severe restriction 
to the models of multiparticle production in proton-proton collisions at high energies. 
This particularly concerns the weighted superposition of two NBDs which is usually 
attributed to the classification of events into soft and semihard with respect to the 
momentum transfer in parton-parton scatterings. 

Figure [1] shows data on MD measured by the ATLAS Collaboration [18] in the 
pseudorapidity window |?7| < 2.5 and the phase-space region p^ > 100 MeV/c, rich > 2 
at the energies y/s = 7 and 0.9 TeV. The distributions are displayed around a maximum 
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(a) (b) 

Figure 1. Data on MD of charged particles in pp collisions (open circles) measured by 
the ATLAS Collaboration [18] in the interval \r]\ < 2.5 for pr > 100 MeV/c and rich > 2 
at (a) ^/s = 7 TeV and (b) ^/s = 900 GeV. The solid and dashed lines represent fits 
using the superposition of two NBDs and the IPPI model with the indicated number 
of components (N), respectively. The black dots correspond to the IPPI model with 
fitted probabilities ai and the restriction a^+i < a^. 



with the insets illustrating their overall shape. The experimental data are depicted 
as open circles and the error bars correspond to the quadratic sum of statistical and 
systematic uncertainties. The solid lines connect the points representing best fits to 
the data using function ([2]) with N=2 and five free parameters. Terms of the fitting 
procedure are explained in more detail in the appendix. The visual contradiction of the 
experimental data at ^/s = 7 TeV (0.9 TeV) with the solid curves for n < 17 {n < 15) 
and the corresponding large values of x^/dof = 128.7/80 {x^/dof = 70.3/46) exclude 
description of the ATLAS data with two NBDs. As can be seen from figure [H the 
maximum of the measured MD is shifted toward higher multiplicities at both energies 
in comparison with the two-component NBD model. 

The dashed lines in figure [U represent the results of fitting data in the framework 
of the IPPI model with the probabilities for i pairs of colliding partons taken in the 
form ai = a*. In that case, the parameter a is given by the maximum number of the 
active parton pairs A^. The IPPI model with three adjustable parameters m, k and A^ 
does not reproduce the ATLAS data well. The minimal value of x^/dof = 2133/82 
obtained for A^ = 6 at = 7 TeV and x^/dof = 117/48 for A^ = 5 at = 900 GeV 
is unsatisfactory. Here we do not include additional conditions in the minimization 
procedure which were applied in [21]. Such conditions impose constraints on the 
parameters m and k of the model and result in even larger values of (1A.1|) . A 
better approximation of the data can be obtained if one allows freedom in the choice of 
the probabilities at while preserving their monotonic decrease ctj+i < aj. This condition 
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follows from the physical picture of the IPPI model characterized by the relations 
Hi = im and ki = ik. With the increase in the number of independent adjustable 
parameters (A^ + 1 in this case), one can get smaller value of /dof = 268/78 for six 
active parton pairs at ^/s = 7 TeV. The minimal value of x^/dof = 91.7/74 at this 
energy is obtained for ten components with the parameters quoted in the appendix. 
The corresponding distribution is depicted by the black dots in figure [Tt^a). The same 
symbols in figure [H^b) represent the best approximation of the data at ^/s = 900 GeV 
with x^/dof = 74.7/46 and with four components taken into account. One can see from 
figure [1] that the distributions within the IPPI model with + 1 adjustable parameters 
show similar discrepancies with the ATLAS data for n < 15 at both energies, -y/i = 7 
TeV and 900 GeV, as it is in the case of superposition of two independent NBDs. 
Somewhat smaller /dof = 91.7/74 for the IPPI model at y/s = 7 TeV is attributed to 
a better description of the data in the tail of the distribution. Regardless of this, neither 
of these models can describe the ATLAS data around the maximum. The experimental 
distributions at maximum are shifted toward higher multiplicities relative to the best- 
fitted parametrizations based on both models. 

The CMS measurements [3] in a similar pseudorapidity window are less critical to 
the mentioned hypotheses. This is illustrated in figure |2] where data on MD of charged 
particles measured by the CMS Collaboration |3] in the interval \r]\ < 2.4 at the energies 
^/s = 7 TeV and 900 GeV are shown. One can see that the CMS data are more or less 
well-described by the two-NBD superposition (solid lines). The corresponding values of 
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Figure 2. Data on MD of charged particles in pp collisions (open circles) measured 
by the CMS Collaboration H in the interval \r]\ < 2.4 at (a) = 7 TeV and (b) 
y/s = 900 GeV. The solid and dashed lines represent fits using the superposition of two 
NBDs and the IPPI model with the indicated number of components (N), respectively. 
The dotted lines correspond to the IPPI model with fitted probabilities and the 
restriction a^+i < a^. 
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X^/dof = 13.4/121 and x^/dof =10.7/62 are at both energies satisfactory. Somewhat 
larger values of x'^/dof = 49.6/124 (x'^/dof = 14.7/64) are obtained within the IPPI 
model (dashed lines) with the probabilities ai = a* and six (four) components at the 
energies y/s = 7 TeV (900 GeV). At closer inspection, however, one can see that the 
dashed line in figure [2](a) deviates from the shape of the experimental distribution for 
8 < n < 35. This can be improved if one allows freedom in the choice of the probabilities 
ai while preserving the condition a^+i < a^. Unlike similar fittings of the ATLAS data, 
an acceptable fit and much better x'^/dof = 6.6/113 can be obtained within the IPPI 
model at y/s = 7 TeV with twelve components and thirteen parameters quoted in the 
appendix. The modification of the parameters has a smaller impact on the fitted 
distribution at lower energy. At ^/s =900 GeV one can get x^/dof = 11.5/63 within 
the IPPI model with three components taken into account. The dotted lines shown in 
figure |5] connect the points which represent best fits to the CMS data using the IPPI 
model with + 1 adjustable parameters. 

The reason for the acceptable values of the least-squares tests are significant 
systematic uncertainties of the CMS data as compared with the experimental data 
presented by the ATLAS Collaboration. Despite the large systematic errors, the values 
of Pn measured with high statistics allow observation of a certain tendency in the CMS 
data. One can notice from figure [2] that open circles representing the mean values of the 
measured probabilities are shifted toward higher multiplicities at the maximum of the 
distribution as compared to the lines depicted in this figure. Moreover, the shift values 
at the maximum are roughly the same as are the corresponding displacements (An ~ 2) 
clearly visible in figure [H The observation supports the consistency of both data in the 
considered sense. On the other hand, the systematic errors of the CMS data influence 
inference from the relatively small values of x^/ dof and can lead to conclusions which do 
not comply with the ATLAS data. As for errors, similar considerations concerning the 
hypothesis of single NBD and the measurements of the UA5 and ALICE Collaborations 
in small pseudorapidity windows were expressed in [8]. 

Motivated by the discrepancies shown above, we study the high statistic data on 
MDs obtained at the LHC in the framework of weighted superposition of three NBDs. 
The aim is to achieve a detailed description of the data including the shape of the 
distributions in the region of maximal values of Pn- The results are compared with two- 
NBD parametrization of the same data. Hereafter we refer to the notion 'three NBD 
superposition' also for the situation if the third NBD is reduced to its Poisson limit with 
^3 = 00. The limiting case of two NBDs plus one Poisson distribution is used in our 
analysis where stated below. 

3.1. Three- component description of ATLAS data 

We have fitted data [18] on MDs measured by the ATLAS Collaboration in the interval 
\ri\ < 2.5 with weighted superposition of three NBDs. The decomposition of the total 
distribution in the most inclusive phase-space region corresponding to the transverse 
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(a) (b) 

Figure 3. MD of charged particles measured by tlie ATLAS Collaboration [18] in the 
pseudorapidity interval \ri\ < 2.5 for px > 100 MeV/c, n > 1 at (a) = 7 TeV and 
(b) y/s = 0.9 TeV. The solid lines represent fits to the data by three-component 
superposition of NBDs. The dash-dot, dash and dash-dot-dot lines show single 
components corresponding to the indicated parameters. 

momentum cut pr > 100 MeV/c and rich > 2 is shown in figures [3]^a) and (b) at the 
energies i/i = 7 TeV and 0.9 TeV, respectively. The experimental data are indicated by 
symbols and the fitted three-component function ([21) is depicted by the solid line. The 
dash-dot, dash and dash-dot-dot lines represent single NBD components of the total 
distribution. The corresponding parameters and values of are quoted in table [H 

One can see from figure [3] and tabled] that the decomposition into three components 
is similar at both energies. The dominant component with the largest probability 



Table 1. The parameters of superposition of three and two NBDs obtained from 
fits to data [18 on MD measured by the ATLAS Collaboration in the pseudorapidity 
window I77I < 2.5 with the cut pr > 100 MeV/c, n > 1 at = 7 and 0.9 TeV. 







yG=7 TeV 






Vs = 0.9 TeV 




i 


ai 


fii 




ai 


Ui 


ki 


1 
2 
3 


n 7Q7+0.062 
-0.089 

0.0801-°?^ 


11 2+°-^ 


J^-'J'J-o.oe 

c 7 +0.9 
-0.8 

24. +1^- 

— 6. 


64+°-20 

^•"^-0.36 

n 25+°-27 

^'•^"3-0. 15 
^^•^^-0.05 


28.5tl^ 

11 0+°-^ 


1.78t°;l 
5.0 t\l 
28. tTi. 



X^/dof = 6.6/(85-8) X^/dof = 5.5/(51-8) 



1 0.412±0.016 11.04±0.16 2.71±0.09 0.774±0.031 12.63±0.44 2.70±0.12 

2 0.588±0.016 39.50±0.58 2.65±0.11 0.226±0.031 32.9 ±1.2 7.53±0.77 



X^/dof = 128.7/(85-5) 



xVdof = 70.3/(51-5) 



A three- component description of multiplicity distributions 



9 




50 100 150 200 20 40 60 80 100 

n n 



(a) (b) 

Figure 4. Normalized residues of MD relative to the weighted superposition of two 
NBDs (Pn^^^) with the parameters listed in table [TJ The symbols correspond to 
data on MD measured by the ATLAS Collaboration [18] in the pseudorapidity interval 
|?7| < 2.5 with the cut pr > 100 MeV/c, n > 1 at (a) = 7 TeV and (b) = 0.9 TeV. 
The lines correspond to fits to the data by three-component superposition of NBDs. 
The insets show the detailed behavior of the residues at low n. 

gives the main contribution aiUi to the total average multiphcity. The other two 
components contribute to the high and low multiplicity region, respectively. The average 
multiplicities of the first and the second component, ni and n2, increase with energy. 
This results in a broadening of the total distribution. The average multiplicity 77,3 ^ 11 of 
the third component is nearly energy independent. Within the errors quoted in tabled! 
the probabilities show weak energy dependence as well. The values of the parameters 
ki increase with decreasing probabilities ctj. The first NBD with the largest probability 
ai is characterized by the smallest parameter ki. The component under the peak of the 
distribution at low multiplicities is narrow with large values of k^. 

We have compared the ATLAS data with the best-fitted weighted superposition 
of two NBDs. The values of the parameters corresponding to the two-component 
hypothesis are stated in table [H The relative residues with respect to the two-NBD 
fits are depicted in figure HI The points correspond to data measured by the ATLAS 
Collaboration in the pseudorapidity interval |?7| < 2.5 for pj- > 100 MeV/c, n > 1 at 
the energy y/s = 7 TeV and y/s = 0.9 TeV. The insets show the detailed structure 
of the residues at low multiplicities. The solid lines are given by the three-component 
description of the data shown in figure [31 One can see from figure H] that the two-NBD 
approximation of the data is unsatisfactory at both energies. The corresponding values 
of x^/dof quoted in tabled] are too large. The high statistic ATLAS data manifest a 
distinct peak around n ~ 11. The description of the peak clearly seen in the residues 
in figure H] was obtained by the third negative binomial component with ~ 11. The 
component is depicted by the dash-dot-dot line in figure |3l 
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Figure 5. MD of charged particles measured by tlie ATLAS Collaboration [18] in the 
pseudorapidity interval \ri\ < 2.5 for px > 500 MeV/c, n > at (a) = 7 TeV and 
(b) y/s = 0.9 TeV. The solid lines represent fits to the data by three-component 
superposition of NBDs. The dash-dot, dash and dash-dot-dot lines show single 
components corresponding to the indicated parameters. 



We have studied the peaky structure of MD using the ATLAS data in the separate 
phase-space region defined by the conditions pt > 500 MeV/ c and n > 0. Figure |5] shows 
the data measured in the window |?7| < 2.5 at ^/s = 7 TeV and 0.9 TeV. The symbols 
and the fines fiave tfie same meaning as in figure [31 The corresponding parameters and 
the respective values of x^/dof are given in table El The parameter takes values 
larger than 100 at ^/s = 0.9 TeV with error comparable with infinity. In such a case we 
have fixed = oo and considered the third NBD as a Poisson distribution. 



Table 2. The parameters of superposition of three and two NBDs obtained from fits 
to data [H] on MD measured by the ATLAS Collaboration in the window I77I < 2.5 
with the cut px > 500 MeV/c, n > at ^/s = 7 and 0.9 TeV. 
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One can see some similarities when comparing the description of the most inclusive 
and the pT-cut biased data shown in figures [3] and [5l respectively. The component with 
the largest probability ai has the smallest parameter ki. Its contribution aifii to the 
total average multiplicity is dominant at both energies. The parameters ki increase as 
the probabilities at decrease. The average multiplicities ni and fi2 increase with energy 
resulting in a broader total distribution at ^/s = 7 TeV. The average multiplicity ^3 ~ 3 
of the third component is nearly energy independent similarly to the most inclusive data 
sample. On the other hand, there are differences between the data shown in figures [3] 
and |5l respectively. With the imposed pt cut, the value of becomes significantly 
smaller relative to the average multiplicities of the first and second component. The 
probability 0:3 of the third component increases with ^/s at the expense of ai in the 
PT-cut data sample. The probability 02 of the second component remains nearly energy 
independent. The increase of 0:3 at ^/s = 7 TeV is accompanied by considerable decrease 
of the parameter k^. Unlike the most inclusive sample, the third component becomes 
more significant (03 ~ 0.28) and wider [k^ ~ 2.5) at -y/i = 7 TeV. These observations 
are supported by the CMS measurements [3] in the window |?7| < 2.4 with the same 
condition pt > 500 MeV/c (see figure fT2] below). 

The normalized residues with respect to the two-NBD parametrization (see table [2]) 
of the ATLAS data measured in the pseudorapidity window \ri\ < 2.5 in the phase- 
space region pt > 500 MeV/c, n > at the energies a/s = 7 TeV and 0.9 TeV are 
depicted in figures [6]^a) and (b), respectively. The data at y/s = 7 TeV show sizeable 
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Figure 6. Normalized residues of MD relative to the weighted superposition of two 
NBDs (P^'^^'^) with the parameters listed in tabled] The points correspond to data 
on MD measured by the ATLAS Collaboration [TS] in the pseudorapidity interval 
I77I < 2.5 for PT > 500 MeV/c, n > at (a) = 7 TeV and (b) Vs = 0.9 TeV. The 
lines correspond to fits to the data by superposition of three NBDs. The insets show 
details of the residues at low n. 
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discrepancies relative to the weighted superposition of two NBDs. The corresponding 
value of x^/dof = 70.4/34 for the two NBD fit is too large. As seen from figure [6](a), 
a clear peak is visible at low multiplicities around n ^ 3. The description of the peak 
was obtained by the third negative binomial component with 77,3 ~ 3 as indicated by 
the dash-dot-dot line in figure [5]^a). The residues at the energy -y/i = 900 GeV shown 
in figure [6](b) are nearly fiat though some indication of a peak around n ~ 3 is visible. 

The experimental analysis [18] of the ATLAS data was complemented by higher 
cuts on multiplicity requiring rich > 20 with pr > 100 MeV/c and rich > 6 with 
Pt > 500 MeV/c, respectively. In the first sample, the high-multiplicity cut does not 
allow to study the third low-multiplicity component within superposition of three NBDs. 
This is evident from figure |3] where the dash-dot-dot line showing the third component 
is negligible for n > 20. The second data sample with the cut rich > 6 allows a check 
of consistency of the three-component description, though leaving all parameters free is 
problematic here as well. Consequently, we have fixed some parameters related to the 
low-multiplicity component in this case. At ^/s = 7 TeV, we have kept the probability 
of the third component at the value = 0.28 obtained for > (see figure [5]^a)). 

Figure [Tl^a) shows data on MD measured by the ATLAS Collaboration at ^/s = 
7 TeV in the pseudorapidity window \ri\ < 2.5 in the phase space region px > 500 MeV/ c 
and n > 5. The data points and the lines are shown in a similar way to figure [51 The 
multiplicity cut n > 5 has an impact on the disappearance of the peak at maximum, 
though its remnant is still visible in the low multiplicity region. This is manifested by 
the increase of P„ as n — )■ 5. The experimental data shown in figure Ul^a.) allow us to 




(a) (b) 

Figure 7. MD of charged particles measured by the ATLAS Collaboration [TS] in 
the pseudorapidity interval \r]\ < 2.5 for pr > 500 MeV/c, n > 5 at (a) ^ = 7 TeV 
and (b) y/s — 0.9 TeV. The solid lines represent fits to the data with the restrictions 
quoted in table |3l The dash-dot, dash and dash-dot-dot lines show single components 
corresponding to the indicated parameters. 
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Table 3. The parameters of superposition of three NBDs obtained from fits to data 
[IH] on MD measured by the ATLAS Collaboration in the window I77I < 2.5 with the 
cut pT > 500 MeV/c, n > 5 at = 7 and 0.9 TeV. 
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determine all free parameters provided 0:3 is fixed. The fitted parameters are quoted in 
figure [Tl^a) and table [31 They have nearly the same values as obtained with the condition 
n > (see figure [5]^a) and table [2]). This result represents a consistency check of the 
description of data samples with different multiplicity cuts. The ATLAS data measured 
in the window \ri\ < 2.5 in the phase-space region px > 500 MeV/c, n > 5 at the energy 
-y/i = 900 GeV are depicted in figure [Tl^b). The multiplicity cut n > 5 does not allow to 
determine any parameter of the low-multiplicity component at this energy. Therefore, 
we have fixed 03 = 0.06, = 3.0 and /C3 = 00 at the values obtained for n > (figure 
[5]^b)) leaving the remaining parameters free. The result of the fit is shown in figure 
[7(b) and table [3l Within the errors indicated, the fitted parameters are more or less 
consistent for both cases with and without the n > 5 cut. 

3.2. Pseudorapidity dependence of the third component and CMS data 

The CMS Collaboration provided the results of systematic measurements [3] of MD of 
charged particles produced in proton-proton collisions in the central interaction region. 
The data were accumulated in five pseudorapidity ranges from \ri\ < 0.5 to \ri\ < 2.4 
at the energies y/s = 900, 2360 and 7000 GeV. The measurements refer to the inelastic 
NSD interactions. The experimental analysis was based on about 132 000, 12 000 
and 442 000 events at = 900, 2360 and 7000 GeV, respectively. The data on 
MD in the central pseudorapidity windows allow us to study the behavior of the fitted 
parameters in dependence on the size of the phase-space regions. We have analyzed 
data recorded by the CMS experiment at y/s = 7 and 0.9 TeV. The reason is that 
the analysis in the framework of the superposition of three NBDs is very sensitive to 
the statistics of data which is sufficiently high at both these energies. The statistics of 
12 000 events are too low to study a three-component description of the CMS data at the 
energy ^/s = 2360 GeV. Though systematic errors of the measurements are relatively 
large especially near the maximum of the distributions, the values of P„, when based 
on high statistics, allow us to extract information on the fitted parameters in most 
cases sufficiently reliably. The large systematic errors of the CMS measurements, as 
compared with the ATLAS data, give relatively small values of even for the two-NBD 
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hypothesis. Following the objectives of this combined analysis, significantly smaller 
can be obtained in larger pseudorapidity windows |?7| < rj^ when superposition of three 
NBDs is considered. The results of the fits in different pseudorapidity intervals are shown 
in table |H Here some guidance was taken from the description of the ATLAS data in 
the window |?7| < 2.5 concerning the values of the parameter listed in table [H Within 
the errors quoted therein, ~ 24 was found to be finite at y/s = 7 TeV. Consequently, 
we have fitted the CMS data using k^ as one of the adjustable parameters in the largest 
intervals \ri\ < 2.4 and |?7| < 2.0 at this energy. The fittings give the corresponding finite 
values of k^ in this region, though uncertainties of their determination from the CMS 
measurements admit description with ^3 = oo. In smaller pseudorapidity windows, 
the mean values of k^ turned out to be even larger than 100 with error intervals of 
infinity. Therefore, we have fixed k^ = 00 for rjc = 1.5, 1.0 and 0.5. As follows from 
table [H the errors of all parameters are larger at smaller energy. Within the errors, the 



Table 4. The parameters of superposition of three NBDs obtained from fits to data [3] 
on MD measured by the CMS Collaboration in the pseudorapidity intervals \ri\ < rjc, 
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Figure 8. MD of charged particles measured by the CMS Cohaboration [7 m the 
pseudorapidity windows \ri\ < ijc, rjc = 2.4, 2.0 and 1.5 at (a) y/s = 7 TcV and (b) 
= 0.9 TeV. The soHd hnes represent fits to the data by the three-component 
superposition of NBDs. The dash-dot, dash and dash-dot-dot hnes show single 
components corresponding to the indicated parameters. 
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(a) (b) 

Figure 9. MD of charged particles measured by the CMS Collaboration [3] in the 
small windows |r;| < 1.0 and \ri\ < 0.5 at (a) = 7 TeV and (b) y/s = 0.9 TeV. The 
solid lines represent fits to the data by three-component superposition of NBDs. The 
dash-dot, dash and dash-dot-dot lines show single components corresponding to the 
indicated parameters. 



value of is consistent with infinity in the window \ri\ < 2.5 for the ATLAS data at 
^/s = 900 GeV. In view of the large systematic uncertainties of the CMS measurements, 
this substantiates the choice ^3 = oo in all pseudorapidity intervals when fitting the 
CMS data at = 900 GeV. 

Figure IS shows the MD of charged particles measured in the pseudorapidity windows 
\t]\ < 2.4, |?7| < 2.0 and |?7| < 1.5 at the energies -y/i = 7 TeV and -y/i = 0.9 TeV. The 
experimental data are depicted by symbols. The solid lines represent fits corresponding 
to weighted superposition of three NBDs. The single components are shown by the dash- 
dot, dash and dash-dot-dot curves, respectively. The zero point Pq was not included in 
the fits. Figure |9] shows the MD measured in the small windows \ri\ < 1.0 and \r]\ < 0.5. 
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The data at the energy y/s = 900 GeV do not allow to extract reliable information on 
superposition of three NBDs in the narrowest window |?7| < 0.5. Here we have fixed 
the probability of the second component at the value 02 = 0.15 obtained in the same 
window at -y/s = 7 TeV and also in the larger window |?7| < 1.0 at y/s = 0.9 TeV. 

One can see from figures [H] and [9] that the three-component description of the CMS 
data in different pseudorapidity ranges manifests analogous trends as observed from the 
analysis of the most inclusive sample of the ATLAS measurements shown in figure [31 
Moreover, values of the fitted parameters in the intervals \ri\ < 2.4 and \ri\ < 2.0 are in a 
sense similar to those obtained from the ATLAS data. The multiplicity component with 
largest probability ai has smallest parameter ki in all windows. Its contributions aifii 
to the total average multiplicities are dominant at both energies. Except for the window 
\r]\ < 0.5 at y/s = 900 GeV one can see that the other two components contribute to the 
high and low multiplicity region, respectively. The values of the parameters ki increase 
with the decreasing probabilities a,. The average multiplicities of the first and the 
second component, ni and n2, increase with energy. The average multiplicity of the 
third component is within errors nearly energy independent for all rjc- The approximate 
energy independence is valid for the probabilities as well. 

We have fitted the CMS data in different pseudorapidity intervals with weighted 
superposition of two NBDs. The obtained values of parameters are listed in table O 
Figure [TO] shows the normalized residues of MD relative to the two-NBD fits at the 



Table 5. The parameters of superposition of two NBDs obtained from fits to data [3] 
on MD measured by the CMS Collaboration in the pseudorapidity windows jr/j < 77c, 
T]c = 2.4, 2.0, 1.5,1.0, 0.5 at = 7 and 0.9 TeV. 
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Figure 10. Normalized residues of MD relative to the weighted superposition of two 
NBDs (P^^^^ ) with the parameters listed in table El The symbols correspond to data 
on MD measured by the CMS Collaboration 3 in the pseudorapidity intervals < -q^, 
ric = 2.4, 2.0 and 1.5 at (a) = 7 TeV and (b) ^/s = 0.9 TeV. The lines correspond 
to fits to the data by three-component superposition of NBDs. The insets show the 
detailed behavior of the residues at low n. 
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Figure 11. Normalized residues of MD relative to the weighted superposition of two 
NBDs (P^^^^) with the parameters listed in table [51 The symbols correspond to data 
on MD measured by the CMS Collaboration I3| in the small windows \ri\ < 1.0 and 
\r]\ < 0.5 at (a) = 7 TeV and (b) y/s = 0.9 TeV. The lines correspond to fits to 
the data by three-component superposition of NBDs. The insets show the detailed 
behavior of the residues at low n. 



energies ^/s = 7 and 0.9 TeV in the windows |?7| < 2.4, |?7| < 2.0 and |?7| < 1.5. The 
residues for small intervals \ri\ < 1.0 and \ri\ < 0.5 are depicted in figure [Til The 
insets represent the detailed structure of the residues at low multiplicity. The solid 
lines correspond to the three-component description of the data shown in figures [8] and 
M Due to large errors, the residual analysis is roughly compatible with hypothesis 
of description of the CMS data by two NBDs in all pseudorapidity intervals at both 
energies. However, in regard to a similar analysis of the ATLAS data one has to look 
at the residues in more detail. Let us note that the error bars shown in figures [10] and 
[TT] form a characteristic envelope around mean values of the residues mainly refiecting 
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systematic uncertainties of the data. The envelope follows the fine structure of the 
residual mean values quite accurately. As seen from figure [TOl the structure indicates 
the systematic emergence of a peak at low multiplicities. The peak in the residues is best 
visible in the large windows |?7| < 2.4 and |?7| < 2.0. A comparison of the corresponding 
insets in figures HI and [TO] shows that the peak has nearly the same properties as observed 
in the ATLAS data. This concerns its position as well as its size. The peaky structure 
at low n diminishes with decreasing width of the pseudorapidity intervals. As seen from 
figure [in the residues in the small windows are flat. At ^/s = 900 GeV, two NBD 
superposition (single NBD) gives a nearly perfect description of the data in the interval 
|?7| < 1.0 (|?7| < 0.5), respectively. Similar holds at the energy ^/s = 7 TeV for \ri\ < 0.5 
where two-NBD superposition approximates the data with high accuracy. 

The CMS Collaboration measured the MD of charged particles in the 
pseudorapidity window \ri\ < 2.4 exploiting the transverse momentum cut pr > 
500 MeV/c [3]. The complementary pt cut represents selection of harder processes 
which has significant impact on the shape of the distribution relative to unbiased data. 
The measurements allow as to study the distribution in the framework of the three- 
component model and compare it with analysis of the ATLAS data obtained under 
the same restriction on the transverse momentum of registered particles. The obtained 
values of the fitted parameters are quoted in table [6] The px-cut data at y/s = 0.9 TeV 
does neither allow to extract reliable information on three-component parametrization 
of MD in the case if the third component is considered in the Poisson form. Therefore, 
we have fixed the probability ^2 = 0.11 at the value found from analysis of ATLAS data 
at the same energy with the same px cut in the interval \ri\ < 2.5 (see table [2]). 

Table 6. The parameters of superposition of three and two NBDs obtained from fits 
to data [3 on MD measured by the CMS Collaboration in the window |?7| < 2.4 with 
the cut PT > 500 MeV/c at = 7 and 0.9 TeV. 
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Figure [T2] shows the MD measured by the CMS Collaboration in the pseudorapidity 
interval \ri\ < 2.4 with the transverse momentum cut pt > 500 MeV/c at the energies 
y/s = 7 and 0.9 TeV. The experimental data are depicted as symbols. The fitted 
superposition of three NBDs is indicated by the solid line. The dash-dot, dash and 
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Figure 12. MD of charged particles in the pseudorapidity window \ri\ < 2.4 measured 
by the CMS Collaboration [3] with the transverse momentum cut px > 500 MeV/c 
at (a) ^/s — 7 TeV and (b) y/s = 0.9 TeV. The solid lines represent fits using three- 
component superposition of NBDs. The dash-dot, dash and dash-dot-dot lines show 
single components corresponding to the indicated parameters. 



dash-dot-dot lines show single negative binomial components of the total distribution. 
A comparison of the results of analysis displayed in figures [12] and |5] leads to similar 
conclusions for both measurements when the cut pt > 500 MeV/c is applied. The 
values of all parameters obtained from fits to the CMS data are compatible with those 
extracted from the ATLAS data. The performed combined analysis shows that the 
average multiplicity of the third component under the peak at low n is approximately 
energy independent. This observation does not depend on the pt cut imposed on the 
data. On the other hand, the condition pt > 500 MeV/c influences the energy behavior 
of the parameter significantly. As seen from figure [12] and table |6l the probability of 
the third component indicates an increasing tendency with energy reflecting growth 
of the peak at low multiplicities at a/s = 7 TeV. The observation is in accord with 
analysis of the ATLAS data with pr > 500 MeV/c (see figure [5] and table [2]). Such 
a trend is nearly negligible in the minimum biased CMS data (figure [H]) and is totally 
absent in the ATLAS data with the low transverse momentum cut pt > 100 MeV/c 
(figure [3]). The obtained value of from the pr-cut CMS data at i/i = 7 TeV does 
not contradict the decreasing tendency of this parameter which is visible in the ATLAS 
data with the same restriction on transverse momentum (figure [5] and table [2]). 

Figure [13] shows the normalized residues of MD with respect to the two-NBD 
parametrization (see table |6|) of the CMS data with the cut pt > 500 MeV/c. Though 
the error bars dominated by systematic uncertainties are large, the residual mean values 
indicate existence of a peaky structure near n ~ 3 — 4 at the energy y/s = 7 TeV. The 
structure is similar to the peak clearly visible in the analysis of the ATLAS data shown 
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Figure 13. Normalized residues of MD relative to the weighted superposition of two 
NBDs (P^^^^ ) with the parameters listed in table |6l The symbols correspond to data 
on MD measured by the CMS Collaboration [3] in the pseudorapidity interval jTyj < 2.4 
with the cut pr > 500 MeV/c at {a) ^ = 7 TeV and (b) Vs = 0.9 TeV. The lines 
correspond to fits to the data by three-component superposition of NBDs. The insets 
show the detailed behavior of the residues at low n. 



in figure EJ^a) . The full line corresponds to description of the CMS data by weighted 
superposition of three NBDs with the third component indicated by the dash-dot-dot 
line in figure [T2ra). The residues at a/s = 900 GeV shown in figure fT3](b) are flat. This 
is in compliance with residual analysis of the ATLAS data for p^ > 500 MeV/c at the 
same energy (figure [6](b)). Both p-p-cut data are here well-parametrized by superposition 
of two NBDs. 

3.3. The third component in the ALICE and LHCb data 

Some structure in the distributions of charged particles in the low multiplicity region 
is seen also in the measurements of two other experiments at the LHC The ALICE 
and LHCb Collaborations presented data which allow to search for complementary 
information on the third component in description of MD in pp collisions at high energies. 

The ALICE Collaboration measured [19] multiplicity of charged particles produced 
in the central pseudorapidity region |?7| < 1 at the energy y/s = 7 TeV. The data 
sample was collected from 240 000 events without normalization to the NSD processes. 
The event class was defined by requiring at least one charged particle in the measured 
pseudorapidity interval. The measurements cover the multiplicity region 1 < n < 70. 
The ALICE data at a/s = 7 TeV are shown as points in figure [T^ a). The solid line 
represents a fitted three-component function ([2]) with the parameters quoted in table [71 
The point Pi was not included in the fit. The parameter ^3 was found to be consistent 
with infinity, having a mean value much larger than 100. Therefore, for better stability 
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Figure 14. (a) MD of charged particles measured by tlie ALICE Collaboration [19] in 
the central window \r]\ < 1.0 at ^/s — 7 TeV. (b) Normalized residues of MD relative to 
the weighted superposition of two NBDs (Pn'^^^) with the parameters listed in table [T] 
The inset shows the behavior of the residues at low n. The solid lines represent fitted 
three-component superposition of NBDs. The dash-dot, dash and dash-dot-dot lines 
show single components corresponding to the indicated parameters. 



of the fit, we set = oo. The single components of the total distribution are depicted 
by the dash-dot, dash and dash-dot-dot lines, respectively. The values of the fitted 
parameters indicate a similar trend as follows from the analysis of the CMS data in the 
pseudorapidity window \ri\ < 1 at the energy y/s = 7 TeV (figure [9]^a)). The average 
multiplicity ni ~ 5 of the third component is approximately the same for both data sets. 
Analysis of the ALICE data gives somewhat smaller values of ni and n2 than obtained 
from the CMS measurements. The parameters ki increase with decreasing probabilities 
Oj. Unlike the distributions measured by the CMS and ATLAS Collaborations, small 
wavy fiuctuations are seen in the ALICE data for n > 25. The fiuctuations are attributed 
to the correlation of errors over a range comparable to multiplicity resolution [19]. The 
wavy behavior does not allow to extract the parameters of the third component of the 
function ([2]) from ALICE data at lower energies. 

We have fitted the ALICE data at = 7 TeV with weighted superposition of 
two NBDs. The obtained parameters are listed in table [3 Figure [TW b) shows relative 
residues of MD with respect to the two-NBD fit to the data. The fitted three-component 
superposition of NBDs is represented by the solid line. The inset depicts the detailed 
behavior of the residues at low n. As follows from figure fMT b). the ALICE measurement 
in the window \ri\ < 1.0 is approximately compatible with the hypothesis of description 
of MD by two NBDs. However, in regard to the residual analysis of the ATLAS data 
complemented by the study of r^c- dependence of the residual structures with the CMS 
data, one needs to look at figure [TW b) in more detail. There are some oscillations of 
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the residues for n > 25 which reflect the fluctuations in data mentioned above. Besides 
this, there is an indication of a little peak around ~ 6. We consider that the peak is of 
the same origin as the residual structure at low multiplicities emerging in the CMS data 
in the window \r]\ < 1 at y/s = 7 TeV (see figure [TlT a)). An indirect support for this 
claim can serve values of the parameters and which are within errors mutually 
comparable for both data sets (see corresponding numbers in tables [Hand [7]). 



Table 7. The parameters of superposition of three and two NBDs obtained from fits 
to the ALICE data 19 in the central region \r]\ < 1.0 and the minimum biased LHCb 
data [25] in the forward region 2.0 < rj < 4.5, all at -^/s = 7 TeV. 
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0.567±0.058 
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0.373±0.087 


24.1 ±2.3 


4.5 ±1.2 



xVdof = 15.7/(65-5-1) 



xVdof = 16.3/(39-5) 



The LHCb Collaboration measured the MD of charged particles at the energy 
^/s = 7 TeV in the forward region [25]. The experimental analysis was based on a 
high statistic sample of 3M events. In this paper we study the data taken with the 
minimum bias trigger defined by requiring at least one reconstructed track in the vertex 
detector. The LHCb data obtained in the forward pseudorapidity region 2.0 < rj < 4.5 
are depicted as points in figure [TSl^a). The solid hne represents the best fit to the data 
by superposition of three NBDs. The single components are shown as dash-dot, dash 
and dash-dot-dot lines, respectively. The parameters of the distribution are stated in 
table [T] The analysis gives a large value of > 100. For the sake of better stability of 
the fit, the parameter was fixed at infinity, ^3 = oo. All parameters ki increase with the 
decreasing probabilities as in the central interaction region. The average multiplicities 
fii of single components match with the values obtained from the analysis of the CMS 
data in the regions with nearly the same widths of the pseudorapidity windows. The 
span of the pseudorapidity interval of the LHCb measurements, Arj = 2.5, corresponds 
to the size of the central windows between |?7| < 1.0 and |?7| < 1.5. All three parameters 
Hi shown in figure [T5](a) are found between the corresponding average multiplicities 
obtained from analysis of the CMS data measured in these windows (see tables H] and 
[7]). The observation complies with the extensive character of the average multiplicity 
for each component separately. 

The shapes of the MD in the forward and central region are similar. One can see 
from figure [TSl^ a) that the minimum biased data measured by the LHCb Collaboration 
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Figure 15. (a) MD of charged particles measured by the LHCb Collaboration [25] 
in the forward region 2.0 < rj < 4.5 at y/s = 7 TeV. (b) Normalized residues of 
MD relative to the weighted superposition of two NBDs (Pn'^^'^) with the parameters 
listed in table [T) The inset shows the behavior of the residues at low n. The solid lines 
represent fitted three-component superposition of NBDs. The dash-dot, dash and 
dash-dot-dot lines show single components corresponding to the indicated parameters. 



in the forward region manifest existence of a narrow maximum and a broad shoulder 
at high muhiphcities. It is worth mentioning that a shght decrease of the probabihty 
a2 accompanies an increasing tendency of the parameter k2- This is connected with a 
receding role of the second component in the tail of the total distribution which becomes 
narrower and suppressed with respect to the dominant one (compare with figure [MT a)). 
On the other hand, the probability as of the third component under the peak at low n 
is nearly the same as in the central region (see figures |9]^a) and[T4]^a)). 

We have fitted the minimum biased LHCb data with weighted superposition of 
two NBDs. The obtained values of the corresponding parameters are listed in table [3 
The normalized residues of MD relative to the two-NBD fit of the data are shown 
in figure ITST b). The inset represents the detailed structure of the residues at low 
multiplicity. The solid lines correspond to the three-component description of the LHCb 
data with the parameters listed in table [71 One can see a certain structure at low 
multiplicities with a local maximum around n ~ 7. The maximum is described by the 
third component of the total MD shown by the dash-dot-dot line in figure ITST a). 

Although the residual structures of MD in the ALICE and LHCb measurements 
are not very significant, the data show signs of a small peak at low multiplicities which 
supports analogous observation in the CMS data in the interval |?7| < 1.0 at y/s = 7 TeV. 
The evolution of the peak with the increasing width of the pseudorapidity window in the 
CMS data and its comparability in the largest windows with the peaky structure clearly 
visible in the ATLAS data represents a sequence of observations indicating emergence 
of a new signal in multiparticle production in hadron collisions at TeV energies. 
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4. Invariant Properties of MDs 

The studies of the MD of charged particles in pp/pp collisions showed a complex structure 
which was attributed by many authors to soft and semihard components in the data. 
This led to the natural conclusion that the soft component prevails at lower collision 
energies while the semihard component should play an increasing role at higher energies. 
The important experimental observation concerns the energy independence of MD of 
the soft events 126] . The energy invariance of the parameters of the soft component 
was confirmed by the recent analysis of the CMS data in the framework of the two-NBD 
model [22]. The results of this study show the stability of the description of the distinct 
peak visible in the region where the soft production processes dominate. The peaky 
structure can be well-described by a separate component within the three-component 
parametrization of data on multiplicities measured at the LHC The component under 
the peak at low n reveals remarkable properties. Its average multiplicity is nearly 
the same at both energies y/s = 7 and 0.9 TeV. For minimum biased data, the third 
component is narrow, well-approximated by the Poisson distribution (or by NBD with 
large k^). Description of the experimental data by superposition of three NBDs manifests 
some invariant properties in the studied kinematic region. These are characterized by 
approximate independence of the corresponding parameters on the collision energy ^/s 
and the size of the pseudorapidity window |?7| < rjc- 

Figure [16] shows the probabilities of single negative binomial components at the 
energies ^/s = 7 and 0.9 TeV as function of rjc- The symbols correspond to the results of 
analysis of the minimum biased CMS data in five windows from rjc = 0.5 to rj^ = 2.4 and 
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(a) 




(b) 



Figure 16. The 7]c dependence of the probabihties ai of single components of weighted 
superposition of three NBDs fitted to the NSD CMS data [5] and the most inclusive 
sample of the ATLAS data [H] at (a) = 7 TeV and (b) = 0.9 TeV. The values 
of as are multiplied by a factor of one third. 
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the most inclusive sample of ATLAS data in the interval \r]\ < 2.5. Despite considerable 
errors, there are clear trends visible in the behavior of the probabilities. All reveal 
signs of independence on the window size at y/s = 7 TeV in the pseudorapidity ranges 
from rjc = 0.5 up to rjc = 2.5. The same property seems to be valid at y/s = 0.9 TeV 
regardless of problems with determination of the parameters in the narrowest window 
at this energy. Moreover, from comparison of figures fTBT a) and (b) one can deduce that 
the probabilities a, are approximately energy independent. 

Figure [T7] shows the average multiplicities Ui of the negative binomial components 
of the total distribution in dependence on the size of the window \ri\ < 77^ at the energies 
^/s = 7 TeV and 0.9 TeV. The values of Ui extracted from fits to the minimum bias 
CMS data are depicted for five consecutive intervals in pseudorapidity from rj^ = 0.5 up 
to rjc = 2.4. The symbols at the point rjc = 2.5 correspond to the parameters obtained 
from analysis of the most inclusive ATLAS data with px > 100 MeV/c and rich > 1- The 
average multiplicities demonstrate an approximate linear increase with the window 
size r]c. Within the errors indicated, the linear dependencies can be extrapolated to 
the point fij = at r/c = 0. These properties refiect the extensive character of the 
average multiplicities for all three components of the total distribution. The energy 
dependences of single are mutually different. While the average multiplicities of the 
first and the second component, Ui and n2, increase with y/s, the average multiplicity 
na of the third component is within errors approximately energy independent for all rjc- 
The last statement is an extrapolation in the pseudorapidity interval |?7| < 0.5 because 
the CMS data at i/i = 0.9 TeV does not allow to extract reliable information on the 
third component in the smallest window. 
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Figure 17. The 77c dependence of the average multiphcities fii of single components 
of weighted superposition of three NBDs fitted to the NSD CMS data [3] and the most 
inclusive sample of the ATLAS data [IH] at (a) = 7 TeV and (b) y/s = 0.9 TeV. 
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Figure 18. The r/c dependence of the inverse values of the parameters ki of single 
components of weig hted superposition of three NBDs fitted to the NSD CMS data [3] 
and the most inclusive sample of the ATLAS data [18] at (a) -y/s = 7 TeV and (b) 
Vs = 0.9 TeV. 



Figure [18] shows the inverse values of the parameters ki of weighted superposition 
of three NBDs at the energies y/s = 7 and 0.9 TeV as a function of r]c. The parameters 
were found from analysis of the same data as in the case of figures [16] and [T7J Within 
the errors indicated, the values of fcj do not depend on r]c at both energies. An exception 
is the window \r]\ < 0.5 at -y/i = 0.9 TeV where ki and /c2 show a tendency to 
decrease. This however should be taken with caution because application of three NBDs 
is problematic in the narrowest window at this energy and the parameters were obtained 
here with fixed value of a2 only. As concerns the energy dependence, the parameter ki 
demonstrates some decrease with -y/s. This means that the dominant component of the 
total distribution reveals broadening as the energy increases. At the same time ^2 and 
^3 are approximately energy independent for the minimum biased data. Therefore, the 
energy dependence of the NBD components in the tail and under the peak of MD is 
given entirely in terms of the respective average multiplicities n2 and ns. Because the 
latter does not depend on energy, the peak at the maximum of the distribution reveals 
invariant properties. The imposed pr cut influences the width of the peak at low n but 
leaves the property of the energy independence of preserved. 

5. Summary 

The MDs of charged particles have been studied in limited intervals of pseudorapidity 
using the high statistic data obtained in proton-proton collisions at the LHC The 
analysis concerns the symmetric pseudorapidity windows with the center ai rj = 0. 
The minimum biased data measured at the LHC in the forward region was analyzed 
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as well. The performed combined analysis of data obtained by the ATLAS, CMS, 
ALICE and LHCb Collaborations shows that besides the shoulder-like structure at high 
multiplicities, MD manifests a distinct peak in the low multiplicity region. The peak is 
best visible in the ATLAS measurements which are based on recording the high statistics 
more than lOM events. The description of the data within two-component superposition 
of NBDs is unsatisfactory. The ATLAS data can be well parametrized by a weighted sum 
of three NBDs. The main features of the three-component description are confirmed 
by the analysis of the CMS data in different pseudorapidity windows. The obtained 
parametrization of MD reveals some invariant properties. There are indications, at 
least at the energy ^/s = 7 TeV, that the parameters ki of single NBDs do not depend 
on the window size. The analysis supports an idea that the pseudorapidity invariance 
concerns the probability of each component as well. The component in the tail of the 
distribution and the component under the peak at low n have nearly energy independent 
shapes for the minimum biased data samples. 

The third component of the total distribution reflects the character of the peaky 
structure of MD at maximum. Its probability is at the level of a few per cent. The 
component is narrow, well-described by the Poisson distribution (or by a NBD with a 
large value of ^3). The corresponding average multiplicity n^, is approximately energy 
independent. The low-multiplicity component does not diminish when a cut on the 
transverse momentum is imposed. Its probability 0:3 increases with the pt cut and 
the respective NBD becomes wider with a smaller value of the parameter k^. The third 
component is best visible in the larger pseudorapidity windows. A sign of the component 
is seen also in the minimum biased data measured by the LHCb Collaboration in 
the forward region. The analysis of the LHC data on MD excludes interpretation of 
this Poisson-like component in terms of hard processes which could lead to 'elbow-like' 
behavior at the tail of the distributions considered in |27] . 

Within the performed analysis and on the basis of the studied material we conclude 
that, besides the shoulder-like structure at high multiplicities, the high statistic data 
on MD measured in pp collisions at the LHC manifest existence of a distinct peak 
visible at low n. The total distribution can be well-described by a weighted sum of 
three NBDs. The parametrization of data reveals remarkable properties concerning the 
energy and pseudorapidity invariance of some characteristics of single components of the 
total distribution. Whether the third multiplicity component corresponds to a specific 
class of events or appearance of the characteristic peak in MD is a manifestation of new 
phenomena in multiparticle production remains an open question. 

We have shown that the ATLAS data on multiplicities measured by recording 
high statistics defy description within the model of IPPIs even if the probabilities of 
single interactions are considered as free parameters. Though the IPPI model takes 
somehow into account correlations of particles via assumed NBDs and their convolutions 
in each binary collision, further correlations between the interacting pairs of partons are 
probably needed to capture the enlarge role of collective effects at LHC energies. A 
special question is whether the observed shift of the data at maximum of MD relative 
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to the distributions shown in figure [T] could be interpreted gnature of a phase 

transition in hadron colhsions at high energies. Investigations in that direction are 
however beyond reach of this paper. Physical interpretation of the obtained results 
represents a challenging problem and requires further experimental and theoretical 
study. 
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Appendix 

The results of our analysis are based on the computation and minimization of fo^^ 
various data sets presented in different forms as concerns the binning in multiplicity 
and error estimation. Let us summarize details of the least-squares criteria and the 
conditions used in the minimization procedure. The value of was computed as 

"12 ( pex pph\2 

E ^ , '"^ (A.l) 

m=mi ^ m 

over the finite range < mi,m2 > of the multiplicity bins. The total number of bins 
considered in experiment is denoted by The value of mi = 1 or 2 depends on 
the condition whether the first bin with the lowest multiplicity is included in the fit or 
not. The bins are defined using the information accessible in the following form: the 
probability in the point and the bin interval (am, 6m). In general, the values 
of not given as integers. In the case where the mth bin consists of a single 

multiplicity n, rim = n and Pf^ = P{n) with P{n) given by (2). For the bins containing 
more multiplicities, such as in the tails of the distributions, we distinguish two cases. If 
the point rim lies in the middle of the bin, = {[bm] + [o^m] + 1)/2 (CMS data), we use 
the average value 

n=[am]+l 

If Um 7^ {[bm] + [dm] + l)/2 as in the ATLAS data, we define 

P^^ = (1 - A)P{[nm]) + AP{[nm] + 1), A = - [um]. (A.3) 
All values of P^ are then renormalized to fulfil the condition 

1712 rn2 

E Pt= E (A.4) 

m=Tni m=m\ 

This ensures the same normalization of the phenomenological and experimental 
distributions over the fitted range of bins < mi,m2 >■ The value of is computed 
according to ( lA.ip from the normalized distributions in each iteration step. The 



experimental errors am are calculated from a quadratic sum of the statistical and 
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systematic uncertainties. For the asymmetric systematic errors we used am = cr^ 

> and am = a" if PP^" < . The results of the analysis with the CERN- 
MINUIT program are given in tables 1-7. The obtained values of jdo] are quoted 
in the separate rows. The number of degrees of freedom, do/, is stated explicitly as 
number of the experimental points minus number of the free parameters and minus the 
first experimental point, if not included into the fit. 

Here we state the parameters of fits of the IPPI model to data discussed in the text. 
We obtained the following values with the condition = a*: /dof =2133/(85-3) for 
m = 15.1, k = 0.926 and = 6 (dashed line in figure [T|^a)); x^/dof =117/(51-3) for 
m = 2.7, k = 3.438 and = 5 (dashed line in figure H^b)); x^/dof =49.6/(127-2-1) for 
m = 7.86, k = 2.103 and = 6 (dashed line in figure EJ^ a)); x^/dof =14.7/(68-3-1) for 
m = 3.39, k = 3.099 and = 4 (dashed line in figure [2](b)). 

Allowing the probabilities to be free parameters with the restriction Oj+i < Oj, 
we obtained: x'^/dof =91.7/(85-11) for m = 3.82, k = 2.896, = 0.4555, ^2 = = 
0.1462, a4 = 0.1069, as = = 0.0477, aj = 0.0331, as = ag = 0.0056 and aio = 0.0055 
(black dots in figure [U^a)); x'^/dof =74.7/(51-5) for m = 4.01, k = 2.868, ai = 0.677, 
aa = 0.158, = 0.157 and = 0.008 (black dots in figure [T](b)); x'^/dof =6.6/(127-13- 
1) for m = 4.184, k = 3.258, ai = 0.4566, = = 0.1750, = 0.0774, ctg = 0.0601, 
ae = 0.0303, a-r = 0.0147, = 0.0088, ag = 0.0006 and aio = an = ai2 = 0.0005 
(dotted line in figure El^a)); x^/dof =11.5/(68-4-1) for m = 4.29, k = 2.833, ai = 0.621, 
a2 = 0.228 and = 0.151 (dotted line in figure [SJ^b)). 
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